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Abstract. I briefly review results obtained within the variational Hamiltonian approach to Yang-Mills theory in Coulomb 
gauge and confront them with recent lattice data. The variational approach is extended to non-Gaussian wave functionals 
including three- and four-gluon kernels in the exponential of the vacuum wave functional and used to calculate the three- 
gluon vertex. A new functional renormalization group flow equation for Hamiltonian Yang-Mills theory in Coulomb gauge 
is solved for the gluon and ghost propagator under the assumption of ghost dominance. The results are compared to those 
obtained in the variational approach. 



INTRODUCTION 

In recent years there have been substantial efforts de- 
voted to non-perturbative studies of continuum Yang- 
Mills theory. Among these are variational solutions of 
the Yang-Mills Schrtidinger equation in Coulomb gauge 
[1-3]. In this approach the energy density is minimized 
using Gaussian-type ansatze for the vacuum wave func- 
tional. In this talk I will first briefly review results ob- 
tained within the Tubingen approach [2] using a modi- 
fied Gaussian wave functional. Then I will present a new 
approach to Hamiltonian quantum field theory, which al- 
lows to use non-Gaussian wave functionals in variational 
calculations [4]. The approach is illustrated by using a 
wave functional with cubic and quartic gluonic terms in 
the exponential and applied to calculate the three-gluon 
vertex. Finally, I will discuss a new functional renormal- 
ization group (FRG) approach to the Hamiltonian formu- 
lation of Yang-Mills theory [5] and present results sup- 
porting the findings of the variational approach [2, 3]. 



VARIATIONAL APPROACH 

Instead of working with gauge invariant wave function- 
als, it is usually more convenient to fix the gauge. A 
particularly convenient choice of gauge is the Coulomb 
gauge V • A = 0, which allows an explicit resolution 



of Gauss' law, resulting in the gauge fixed Yang-Mills 
Hamiltonian [6] 

H ™ = \J d D x(j- 1 [A}IV[A}n + B 2 )+H c , (1) 

Hc= 2" / dD ^ x ') J ' 1 { A ]p a ( x ) J [A]F ab (x,y)p b (x), 

where 11" (x) = 8/i8A a (x) is the canonical momentum 
(electric field) operator and 

J[A] = Det(-DV) (2) 

is the Faddeev-Popov determinant. Furthermore 

p a ( x ) = -f bc A b W (3) 

is the color charge of the gluons and 

F» b (x,y) = (^al(-DV)- 1 (-V 2 )(-DV)- 1 b,fe) (4) 

is the so-called Coulomb kernel. In the presence of mat- 
ter fields with color charge density p" n {x), the gluon 
charge p a (x) in the Coulomb term is replaced by the 
total charge p a (x) + p^(x) and the vacuum expectation 
value of F ab (x,y) acquires the meaning of the static non- 
Abelian Coulomb potential, see eq. (16) below. 

The gauge fixed Hamiltonian eq. (1) is highly non- 
local due to the Coulomb kernel F ab (x,y), eq. (4), and 
due to the Faddeev-Popov determinant J[A], eq. (2). In 
addition, the latter also occurs in the functional integra- 
tion measure of the scalar product of Coulomb gauge 
wave functionals 



1 Invited talk given by H. Reinhardt at "T(r)opical QCD 2010". 
September 26-October 1, 2010, Cairns, Australia. 



(Y l \0\y 2 ) = jDAJ\A}\rflA]Oy 2 \A}. (5) 



While the elimination of unphysical degrees of freedom 
via gauge fixing is often beneficial (and sometimes un- 
avoidable) in practical calculations, the price to pay is 
the increased complexity of the gauge fixed Hamiltonian 
eq. (1). The non-trivial Faddeev-Popov determinant re- 
flects the intrinsically non-linear structure of the space of 
gauge orbits and dominates the IR behavior of the theory. 
Once Coulomb gauge is implemented, any functional of 
the (transverse) gauge field is a physical state. 

Treating the Hamiltonian (1) in the ordinary Rayleigh- 
Schrodinger perturbation theory one finds to leading or- 
der the usual one-loop j3 function with j3o = —^-N c [7]. 
Our main interest lies, however, in the IR sector of the 
theory, which requires a non-perturbative treatment. 

In ref. [2] the vacuum energy density was minimized 
using the following ansatz for the vacuum functional 



: exp 



d*k 



(0(*)A?(*)A?(-*) 



(6) 

For this wave functional the static gluon propagator is 
given by 1 



(A?(-k)A b j (q)) = S ab t ij (k)(2nfS(k-q) 



2(0 (Jk) 



(7) 



so that (0(k) represents the (quasi-) gluon energy. 

To the order considered in ref. [2], i.e., up to two 
loops in the energy (which corresponds to one loop in the 
associated Dyson-Schwinger equations), it was shown 
[8] that the Faddeev-Popov determinant, eq. (2), can be 
represented as 



J [A] = exp 
where 



- jd\d i yAl(x)xtf(x 1 y)A b ] (y) 



(8) 



1 / 8 2 In J[A] 

2\SAf(x)SA b j(y) 



= 8 a » ti j(x)x(x,y) 



(9) 

is the ghost loop and represents a measure of the curva- 
ture [2] of the Coulomb gauge fixed configuration space. 
Expressing the ghost Green's function as 



((-DV)- 1 } 



S(-A) ' 



(10) 



where d(k) is the ghost form factor, its Dyson-Schwin- 
ger equation reads to one-loop order 



d- l (k)=g- 1 -I d (k) 1 



(11) 



1 Here and in the following, fy(fc) = Sjj —kjkj/k 2 denotes the trans- 
verse projector in momentum space. 



I d (k) = 



Nr 



d q 



ft\2 



i-(M) 



d(q-k) 

<»(«) (q- k ) 2 ' 



(12) 



The ghost form factor d{k) measures the deviation from 
the QED case, where d(k) = 1 . Furthermore, d~ l (k) rep- 
resents the dielectric function of the Yang-Mills vacuum, 
and the so-called horizon condition 



d-\0) = 



(13) 



ensures that the Yang-Mills vacuum is a dual supercon- 
ductor [9]. Minimization of the vacuum energy (H Y m) 
with respect to the kernel co(k) leads to the gap equation, 
which after renormalization reads [2], [10] 

co 2 (k)=k 2 + X 2 (k) +c 2 + A/P) (*) + 

+ 2 X (k){AI (l Hk)+c l }, (14) 

where 

A/W(*) =/(")(*) -/W(0), (o(Jfc) = co(k)-x(k), 



and 



V(k)^g 2 (F) = 



W n (q)-m n (k) 
(0(q) 

d 2 (k) 



(15) 



(16) 



is the non-Abelian Coulomb potential, where the approx- 
imation 

((-DVrM-V^-DV)- 1 ) 

»{(-DV)- l )-{-V 2 ).{{-DV)- 1 ) (17) 

has been used. In the gap equation (14), c\ and C2 are 
(finite) renormalization constants. For the critical solu- 
tion, where one imposes the horizon condition (13), both 
co(k) and %(k) are infrared divergent, which implies that 
the transverse gluon propagator vanishes at k — > 0, while 



(5(0) =lim(fl>(*) -*(*)) =ci. 

k^fQ 



(18) 



A perimeter law of the 't Hooft loop requires c\ = 
and this value is also favored by the variational principle 
[10]. Furthermore, for ci =0 the IR limit of the wave 
functional (6) is 



W \A} = ^I\y(k), 

k 



v(* = o) = i, 



(19) 



which is the exact vacuum wave functional in 1 + 1 
dimensions [11]. The renormalization parameter C2, on 
the other hand, has no influence on the IR or U V behavior 
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FIGURE 1. Gluon propagator obtained with a Gaussian 
(dashed line) and a non-Gaussian functional (straight line), 
compared to the lattice data from Ref. [12]. 



of the solutions of the gap equation (14). Only the mid- 
momentum regime of (0(k) is weakly dependent on 
[2]. Since we are mainly interested in the IR properties 
we will put C2 = 0. 

The coupled integral equations (11) and (14) can be 
solved analytically in the IR (for the critical solution 
satisfying the horizon condition (13)) using power law 
ansatze [2, 13] 

fi)(Jt) ~ k~ a , d{k)~k-P. (20) 

Under the assumption of a trivial scaling of the ghost- 
gluon vertex one finds that in D space dimensions the IR 
exponents satisfy the sum rule 

a = 2j3+2-£>. (21) 

This sum rule ensures that % (Jfc) has the same IR behavior 
as 0)(k), cf. eq. (18). Furthermore, in D = 3 one finds two 
solutions [13] 

j3~0.8 and J3 = 1. (22) 

The latter gives rise to a Coulomb potential (16) which 
is strictly linear at large distances. In D — 2 only one 
solution J3 = 5 is found, implying a = 5 and a Coulomb 
potential (16) rising as V(r) <~ r 4 / 5 [14]. Furthermore 
in D = 2 only the critical solution satisfying (13) exists 
[14]. On the lattice, on the other hand, one finds a = 1 
in D = 3 and a linearly rising potential [15]. In the UV- 
regime the coupled equations (11) and (14) can be solved 
perturbatively and one finds at large momenta k 3> 1 [2] 

co(k) ~ k, d{k) ~ l/Vlog(£). (23) 

The full numerical solutions of the above equations were 
given, for D = 3 space dimensions, in refs. [2] and [3], 
where the critical solutions j3 = 0.8 and J3 = 1, respec- 
tively, were found, and for D = 2 in ref. [14]. One finds 
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FIGURE 2. The numerical solution for the static non- 
Abelian Coulomb potential [3]. 

an inverse gluon propagator which in the UV behaves 
like the photon energy but diverges in the IR, signal- 
ing confinement in agreement with the IR analysis (20), 
(21), (22). Figures 1 and 2 show the resulting gluon en- 
ergy co(k) and non-Abelian Coulomb potential (16) for 
the solution J3 = 1. Their IR behavior is in agreement 
with the results of the IR analysis. The obtained prop- 
agator also compares favorably with the available lat- 
tice data [12]. There are, however, deviations in the mid- 
momentum regime (and minor ones in the UV) which 
can be attributed to the missing gluon loop, which es- 
capes the Gaussian wave functionals. These deviations 
are presumably irrelevant for the confinement properties, 
which are dominated by the ghost loop (which is fully in- 
cluded under the Gaussian ansatz), but are believed to be 
important for a correct description of spontaneous break- 
ing of chiral symmetry [16]. 

NON-GAUSSIAN WAVE FUNCTIONALS 

In ref. [4] the variational approach to Yang-Mills theory 
in Coulomb gauge was extended to non-Gaussian wave 
functionals of the form 

|y/[A]| 2 =exp(-S[A]), (24) 

S[A] = J coA 2 + 1 J 73A 3 + 1 J 74A 4 , (25) 

where (0, 73, 74 are variational kernels. 

Representing the Faddeev-Popov determinant by a 
functional integral over ghost fields from the identity 

= J DA -^( J {A}e- S[A] K{A}), (26) 

with K[A] an arbitrary functional of the gauge field, one 
can derive, in the standard fashion, DSEs for the various 
gluon and ghost Green functions. These DSEs are the 
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FIGURE 3. Truncated DSE for the three-gluon vertex, under 
the assumption of ghost dominance. 



usual DSEs of Landau gauge Yang-Mills theory, how- 
ever, in D = 3 dimensions and with the bare vertices of 
the usual Yang-Mills action replaced by the variational 
kernels (0, 73, 74. It should be stressed that these Hamilto- 
nian DSEs are not equations of motion in the usual sense, 
but rather relations between the Green functions and the 
so far undetermined variational kernels. By using these 
DSEs, the energy density can be written as a functional 
of the variational kernels, 



(H YM ) =£[©,73,74] . 



(27) 



By using a skeleton expansion, the vacuum energy can be 
written at the desired order of loops. Confining ourselves 
to two loops, the variation of the vacuum energy Eq. (27) 
with respect to the kernel 73 fixes the latter to 



y?£(p,q,k) = 2igf b < 

SjjiP ~ qh + S jk (q - k)j + Skjjk - p)j 
Q.{p)+Q.{q) + a{k) 



(28) 



Equation (28) is reminiscent of the lowest-order pertur- 
bative result [17], with the perturbative gluon energy \p\ 
replaced by the non-perturbative one £l(p). 

With this result, variation of (H YM ) with respect to (0 
yields the gap equation (14), which now, however, con- 
tains on the r.h.s. in addition the gluon loop I g (k) [4], 
which escapes the Gaussian wave functional (6). The 
presence of the gluon loop I g (k) in the gap equation mod- 
ifies the UV behavior and allows us to extract, from the 
non-renormalization of the ghost-gluon vertex, the cor- 
rect first coefficient of the j3 function. In order to esti- 
mate the size of the gluon-loop contribution to the gluon 
propagator, we use the gluon and ghost propagators ob- 
tained with a Gaussian wave functional [3] to calculate 
the gluon loop. The result is shown in Fig. 1, together 
with lattice data from Ref. [12]. The agreement between 
the continuum and the lattice results is improved in the 
mid-momentum regime by the inclusion of the gluon 
loop, i.e., the three-gluon vertex, as observed also in Lan- 
dau gauge [18]. 

The truncated DSE for the three-gluon vertex T3 under 
the assumption of ghost dominance is represented dia- 
grammatically in Fig. 3. Possible tensor decompositions 
of the three-gluon vertex are given in Ref. [19]. For sake 
of illustration, we confine ourselves to the form factor 
corresponding to the tensor structure of the bare three- 
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FIGURE 4. Form factor fy± of the three-gluon vertex for 
orthogonal momenta and comparison to lattice data for the d = 
3 Landau-gauge vertex [20]. The momentum scale is arbitrary 
and has been adjusted to make the sign change occur at the 
same point. The lattice data are shown by courtesy of A. Maas. 



gluon vertex 



where Y^ is the perturbative vertex, given by Eq. (28) 
with Cl(p) replaced by \p\. Furthermore, we consider a 
particular kinematic configuration, where two momenta 
have the same magnitude 



Y^Y 



(0) 



r (0) 

1 3 



AO)- 
1 3 



(29) 



-2 _2 „2 
P\ = Pl=P , 



Pi -P2 = CP 



(30) 



To evaluate the form factor /3a(p 2 ,c), we use the ghost 
and gluon propagators obtained with a Gaussian wave 
functional [3] as input. The IR analysis of the equation 
for fi A (p 2 ,c) [Eq. (29)] performed in Ref. [13] shows 
that this form factor should have a power law in the IR, 
with an exponent three times the one of the ghost dress- 
ing function; this is confirmed by our numerical solution 
[4]. The result for the scalar form factor f^A for orthogo- 
nal momenta, /(/? 2 ,0), is shown in Fig. 4, together with 
lattice results for d = 3 Landau gauge Yang-Mills the- 
ory. Our result and the lattice data compare favorably in 
the low-momentum regime. In particular, in both stud- 
ies, the sign change of the form factor occurs roughly at 
the same momentum where the gluon propagator has its 
maximum. (The scale in Fig. 4 is arbitrary.) 



HAMILTONIAN FLOW 

The advantage of the variational approach to the Hamil- 
tonian formulation is its close connection to physics. The 
price to pay is the apparent loss of manifest renormal- 
ization group invariance. Renormalization group invari- 
ance is naturally built-in in the functional renormaliza- 



tion group approach to the Hamiltonian formulation of 
Yang-Mills theory proposed in ref. [5], 

In the FRG approach the quantum theory of a field <p 
is infrared regulated by adding the regulator term 

AS k [(p] = ^(p-R k -(p = ^J(pR k (p (31) 

to the classical action, which in the Hamiltonian ap- 
proach is given by the logarithm of the wave functional 
(24). The regulator function R k (p) is an effective mo- 
mentum dependent mass with the properties 



Km R k (p) > 0, 

p/k->0 



lim R k (p) = 0, 



(32) 



which ensures that R k (p) suppresses propagation of 
modes with p < k while those with p > k are unaffected 
and the full theory at hand is recovered as the cut-off 
scale k is pushed to zero. Wetterich's flow equation for 
the effective action T k [(j)] of a field (j) is given by 



kd k r k [<j>] = -Tr 



2 tf[4>]+R, 



kd k R k , (33) 



4n) 
kX 



M = 



8"T k [<j>] 



(34) 



50! . . . 8<p n 

are the one-particle irreducible n-point functions (proper 
vertices). The generic structure of the flow equation (33) 
is independent of the details of the underlying theory, but 
is a mere consequence of the form of the regulator term 
(31), being quadratic in the field. By taking functional 
derivatives of Eq. (33) one obtains the flow equations 
for the (inverse) propagators and proper vertices. The 
general flow equation (33) still holds for the Hamiltonian 
formulation of Yang-Mills theory provided that (j) is 
interpreted as the superfield (j) — (A,c,c). 

The FRG flow equations embody an infinite tower of 
coupled equations for the flow of the propagators and the 
proper vertices. These equations have to be truncated to 
get a closed system. We shall use the following trun- 
cation: we only keep the gluon and ghost propagators, 
which we write in the form 

„2 



.(2) 
Lcc 



dkiP) 



(35) 



(3) 

In addition, we keep the ghost-gluon vertex Ti { - , which 
we assume to be bare, i.e., we do not solve its FRG 
flow equation. Moreover, we shall assume infrared ghost 
dominance and discard gluon loops. Then the resulting 
flow equations of the ghost and gluon propagator reduce 
to the ones shown in Figs. 5 and 6. 

These flow equations are solved numerically using the 
regulators 



R A ,k{p) = 2pr k (p), R c>k (p) = p 2 r k (p), (36) 
n 2 k 2 



rk(p) = exp 



(37) 
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FIGURE 5. Truncated flow equation of the gluon propagator. 
The bare vertices at k = A are symbolized by small dots and the 
regulator insertion kd k R by a crossed square. 
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FIGURE 6. Truncated flow equation of the ghost propagator 



and the perturbative initial conditions at the large mo- 
mentum scale k = A, 



d\(p) = d\ = const. 



®a{p) = p + a. 



(38) 



With these initial conditions, the flow equations for the 
ghost and gluon propagators are solved under the con- 
straint of infrared scaling for the ghost form factor. The 
resulting full flow of the ghost dressing function is shown 
in Fig. 7. 

As the IR cut-off momentum k is decreased, the ghost 
form factor d k (p) (constant at k = A) builds up infrared 
strength and the final solution at k = k m i„ is shown in Fig. 
9 together with the one for the gluon energy (O kmjn (p) in 
Fig. 8. It is seen that the IR exponents, i.e., the slopes of 
the curves d kmin (p), (O k ,„ in (p) do not change as the min- 
imal cut-off k m j„ is lowered. Let us stress that we have 
assumed infrared scaling of the ghost form factor but 
not the horizon condition d^Jp = 0) = 0. The latter 
was obtained from the integration of the flow equation 
but not put in by hand (the same is also true for the in- 
frared analysis of the Dyson-Schwinger equations fol- 
lowing from the variational Hamiltonian approach, i.e., 
assuming scaling the DSEs yield the horizon condition). 

The infrared exponents extracted from the numerical 
solutions of the flow equations are a = 0.28 and j3 = 
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FIGURE 7. Flow d k (p) of the ghost form factor. 
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FIGURE 8. Inverse gluon propagator ft) at three minimal cut- 
off values k m i n . 
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FIGURE 9. Inverse ghost form factor d at three minimal cut- 
off values k m i n . 

0.64. They satisfy the sum rule found in [13] resulting 
from the DSE obtained in the variational approach but 
are smaller than the ones of the DSE. Moreover, the 
present approach allows to prove the uniqueness of the 
sum rule (21) [5], analogously to the proof in Landau 
gauge [21]. 

Replacing the propagators with running cut-off mo- 
mentum scale k under the loop integrals of the flow equa- 
tion by the propagators of the full theory, 

d k {p) -> d k=0 (p), co k (p) -> co k=0 (p), (39) 

amounts to taking into account the tadpole diagrams [5]. 
Then the flow equations can be analytically integrated 
and turn into the DSEs obtained in the variational ap- 
proach [2], with explicit UV regularization by subtrac- 
tion. 2 This establishes the connection between these two 



2 Instead of the complete right-hand side of Eq. (14) one really obtains 
only the IR-dominant contribution % 2 - However, it turns out [5] that the 
numerical solution of the equations is the same as in Ref. [2] over the 
whole momentum range. 



approaches and highlights the inclusion of a consistent 
UV renormalization procedure in the present approach. 
The above results encourage further studies, which in- 
cludes the flow of the potential between static color 
sources as well as dynamic quarks. 
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